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Abstract 

We derive the renormalization group equations (RGE) for the flavour coupling 
matrices of the effective dimension-five operators which yield Major ana neutrino 
masses in the multi-Higgs-doublet Standard Model; in particular, we consider the 
case where two different scalar doublets occur in those operators. We also write 
down the RGE for the scalar-potential quartic couplings and for the Yukawa cou- 
plings of that model, in the absence of quarks. As an application of the RGE, we 
consider two models which, based on a \x-t interchange symmetry, predict maxi- 
mal atmospheric neutrino mixing, together with U e z = 0, at the seesaw scale. We 
estimate the change of those predictions due to the evolution of the coupling ma- 
trices of the effective mass operators from the seesaw scale down to the electroweak 
scale. We derive an upper bound on that change, thereby finding that the radiative 
corrections to those predictions are in general negligible. 
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1 Introduction 



The Standard Model (SM) in the strict sense, i.e. without right-handed neutrino singlets, 
forbids neutrino masses. However, it was noticed a long time ago pQ that, if one allows 
for lepton- number nonconservation, then one can construct, with the SM multiplets, op- 
erators of dimension higher than four which give Majorana masses to the neutrinos. The 
lowest-dimensional such operators have dimension five and contain two left-handed lepton 
doublets and two Higgs doublets; those operators can be thought of as arising from the 
seesaw mechanism |2] after one has integrated out the right-handed neutrino singlets. 1 
Under the assumption that the SM is valid up to the seesaw scale rriR, the renormaliza- 
tion group evolution of the dimension-five neutrino mass operators from rriR down to the 
electroweak scale, as represented for instance by the Z° mass mz, can be determined; the 
evolution equations have been computed in the SM and in its minimal supersymmetric 
extension jH El El E] (for a review see (HI). This is an important issue in view of testing 
mechanisms and symmetries for explaining the neutrino masses and the lepton mixing 
angles, since such mechanisms and symmetries are usually operative or imposed at the 
seesaw scale, while the measurements are effected at the electroweak scale. (For the ex- 
perimental and theoretical status of neutrino masses and lepton mixing see, for instance, 
[Oj and [inj, respectively.) 

In this paper we extend the existing results for the SM renormalization group equations 
(RGE) to the case of an arbitrary number of Higgs doublets. In particular, we focus 
on dimension-five neutrino mass operators which contain two different Higgs doublets; 
indeed, to our knowledge, that case has not yet been treated in the literature. The reason 
to consider the multi-Higgs-doublet SM is that, within that framework, several models 
have been produced in recent years which predict, for instance, lepton mixing angles 
013 = and 6 23 = vr/4 [TU d, or #13 = alone CHI El , or 6 23 = vr/4 and 5 = vr/2 H^. 2 
(See [HIGH], for instance, for the definition of the lepton mixing angles.) Those predictions 
usually hold at the seesaw scale and, in order to compare them with experiment, one needs 
to know the corresponding corrections at the electroweak scale. 

In Sect. |21 we display the Lagrangian of the multi-Higgs-doublet SM, without quarks 
but with dimension-five neutrino mass operators, and present the RGE for the couplings of 
that Lagrangian. In Sect.EJwe discuss the specific RGE for the models, referred to above, 
which predict 6*13 = and #23 = 7r/4 at the seesaw scale. In Sect.EJwe show explicitly how 
those sessaw-scale predictions arise, and how they may be changed by the renormalization 
group evolution. In Sect. El we derive an upper bound on the effect of that evolution. A 
short summary of our main results is provided in Sect. El An appendix contains some 
details of the calculation of the beta functions for the neutrino mass operators. 

1 The effect of the dimension-six operators which also arise from the seesaw mechanism has been 
studied in 0. 

2 The predictions #23 = ir/4 and 6 = tt/2 have first been obtained in a supersymmetric extension of 
the SM 118 . 
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2 General case 



2.1 The model 

We consider the SM with n# Higgs doublets 4>i = 1> 2, . . . , nn) with weak hypercharge 
1/2. The SU{2) gauge coupling constant is denoted g while the U(l) gauge coupling 
constant (with the above normalization for the weak hypercharge) is denoted g' . The 
scalar potential V has the form 

V = quadratic terms + ^yW {$,4>j) (filfii) > (!) 

i,j,k,l=l 

where the dimensionless couplings Ayjy satisfy 

Ayfei = Afcjy = X* ilk . (2) 
The lepton Yukawa Lagrangian Cye is 

£vt = - E (M^l + DjXUiiR) , (3) 

i=l 

where £)l denotes the left-handed lepton doublets and £r the right-handed charged-lepton 
singlets. We have defined the dimensionless flavour coupling matrices Yj in the same way 
as [HIE]- Note that in this paper we do not use the summation convention. 
The effective dimension-five neutrino mass operators are defined as 

0«= E E {Dl aa ^C^D hpc ) {e a % b ) , (4) 

a,/3=e,fi,r a,b,c,d=l 

where, contrary to what we had done in (0) and (J3J), we have made explicit both the 
flavour and gauge-S , f/(2) indices, and the summations thereover. In Pj). C is the Dirac- 
Pauli charge-conjugation matrix; a and /3 are flavour indices; a, b, c, and d are SU(2) 
indices; and e is the antisymmetric 2x2 matrix, with e 12 = 1. The flavour coupling 
matrices k^' in (jlj) have dimension —1 and satisfy 

«S ) = *S£?, i-e- « W) = « tj ' i)T . (5) 

2.2 The RGE 

The RGE are first-order differential equations which give the evolution of the couplings 
of a model relative to t = In /i, where /x is the mass parameter used in the regularization 
of ultraviolet-divergent integrals; the basic equation is 

d A k 1 it 



(2ir) (k 2 — m 2 ) (k+p) — m 



where p is a typical momentum and m, m! are typical masses appearing in the loop 
integral; the dots represent either /i-independent terms or terms which disappear in the 
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limit that regularizes the integral. We have computed, at the one-loop level, the RGE for 
the model outlined in the previous subsection. For the RGE of the coupling matrices of 
the effective mass operators we have found 

167T 2 S±- = -3A fe) + 4 E W (fc ° + E [Tki* m + + k™p + p t k^ 
dt fej=i fc=i 



n H 

"t 



+2 e {/c^y/n - 



k=l 



K (ik) + K («) 



A: 



+Y^Y*k^ - YfY* + k^] } , (7) 



where 



T l3 : = tr^y/ 



i in 

p : = sE^n. (9) 

z fe=i 

The third line of (jJJ) is obtained from the second line through the interchange % j 
together with transposition, in agreement with ((HJ). Our result coincides, when i = j, 
with the result given in [HI El; it generalizes that result for the case % ^ j. Note that, for 
the sake of simplicity, in the present paper we dismiss quarks; in general, one would have 
to add to (JHJ) analogous trace terms featuring the Yukawa-coupling matrices of the Higgs 
doublets (pi and <pj to the up and down quarks, multiplied by a colour factor 3 — see, for 
instance, 

The terms in the second and third lines of (J7J) arise from diagrams like the one in 
figure H We dwell on the explicit derivation of those terms in the appendix. 

In order to solve the RGE for the effective neutrino mass operators one also needs 
the RGE for the other couplings occurring in (|7|i . The general RGE for an arbitrary 
renormalizable gauge field theory have been derived in ^3 ^] at the one- and two-loop 
levels, respectively. It is convenient to have the results of jTTj specialized to the case of 
the multi-Higgs-doublet SM. We have found that 

dA ■ UH 

167T — 4 ^ ] (2 ^ijmn^nmkl ^ijmn^kmnl ^imnj^mnkl 

at m,n=X 

~\~^imkn^mjnl ~\~ ^mjkn^imnl ^ijkl 

H ~ Oijdki H — {2dud kj - dijdki) 

<5 4 

^ ] (T m j y^imkl Tml^ijkm Tim^mjkl l^km^ijml) 
m=l 



-2tT{YtfY k Y?), (10) 

16vr 2 ^ = E (T*Y k + Y k YiY + \ Y^ k Y k ) - ^ ^ Y t . (11) 

It is well known that the RGE for g and g' are 

2 d£f / 22 AN n H \ 3 



16tt 2 -^ = + + — )g, (12) 

dt V 3 3 6 1 ' 
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Figure 1: A typical vertex correction in the renormalization of the operator Oij. The 
relevant Yukawa-coupling matrices are indicated. 



_ 2 dg' (20N , n H \ _ /3 



167r 2 ^- = — + (13) 
dt V 9 6 1 



where N = 3 is the number of fermion families. 



3 Application of the RGE to two particular models 

3.1 The 1 9 and Da models 



We now apply the general RGE derived in Sect. El to the so-called Z 2 [TjQj and D4 [12] 
models — for a review see [T9*j . Those models predict 

^13 = 0, 

^23 = vr/4 li4j 

at the seesaw scale and are, in what regards the practical application of the RGE, identical. 
They both have three Higgs doublets 0i, 2 , and 3 . Below the seesaw scale the structure 
of both models is dictated by the symmetries 

ry(aux) . 



J 2 



eR -> -e R , 0i -> -0i, (15) 

: £>L M <"> Dhr, /iR «-> T R , 03 -> -0 3 . (16) 

These two symmetries hold in between the seesaw (high) scale m# and the electroweak 
(low) scale m z . Indeed, they are broken only spontaneously, by the vacuum expectation 
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values (VEVs) of 05 and 0°, respectively, at the low scale. Because of the symmetries 
in (JT5|l and (fTBjl. the Higgs potential is given by 



V 



2 

quadratic terms + A, [4>\<i>?j 
i=i 

+A 4 (<M0i) (4fo) + As (0 f i0i) + Ae (40 2 ) (40; 

+A 7 (0l0 2 ) (0 2 0x) + A 8 (4<i>i) + ^9 (0^ 3 ) (0 3 2 ) 



+ 



A10 (0i0 2 ) 2 + An (0i03) 2 + Ai 2 (0^3) 2 + H.c. 



(17) 



where A10, An, and A12 are the only non-real quartic couplings. Comparing (fTTjl with ((TJ) 
and (j2J), we arrive at the identifications \uu = Aj (for i = 1,2,3), A1122 = A2211 = A4/2, 
A1221 — A2112 = A7/2, A1212 = A2121 = A10, and so on. 

The 7j2 and -D4 models have three other symmetries, the family-lepton-number sym- 
metries L a , which are broken at the seesaw scale — softly in the 7L 2 model, spontaneously 
in the -D 4 model. Because of the symmetries in (|15|) and (|16|). and also because of the 
symmetries L a — which remain valid for the quartic couplings of the light fields below the 
seesaw scale — the lepton Yukawa Lagrangian is 



Y£ 



Le 



1/4 



(D Lfl fi R + £> Lt Tr) 02 - 2/5 (Dl^R ~ -DltT-r) 03 + H.c. (18) 



(The coupling constants y\ y2 occur in the Yukawa interactions of the right-handed neutrino 
singlets and are thus of no concern here.) Comparing (fTBj) with ©, we see that the 
Yukawa-coupling matrices are 



Hence, from 



Y-i 


= diag 


(vl 0, 


0), 


Y 2 


= diag 


(0, yt 


2/1), 


Y 3 


= diag 


(0, yt 


-yl) 




Tn = 


1 |2 

m\ , 






T 22 = 


2M 2 


i 




T33 = 


2N 2 


5 



(19) 



(20) 
(21) 
(22) 

and the T^- with i ^ j vanish, a fact which simplifies considerably the RGE in this 
particular case. 

As emphasized before, the symmetries (II 5|) and (fTT)|) are broken only at the electroweak 
scale. The validity of the symmetry Ti^ — which changes the sign of 0i but does not 
affect the lepton doublets D L — has an important consequence: the operators Oi 2 , 2 ±, 
(9i3, and 03i are altogether absent. The symmetry Z 2 tr \ on the other hand, changes the 
sign of 03 simultaneously with the interchange of with D^ T . This implies that the 
coupling matrices (i = 1, 2, 3) must be of the form 



(23) 
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while the matrices and = are of the form 

/ p -p \ 

q s r 

\ —q — r — s J 



(24) 



3.2 The RGE for the 1 2 and D A models 

The renormalization group equations for the Yukawa couplings of the Z 2 and _D 4 models 
are 

2/3, 



2 dy 3 (5^ |2 9^ + 150' 



167T 2 ^ = l-\y 3 
dt V 2 |y3 



167T 



dy4 
dt 



I J/4 



4 

I |2 
2/5 - 



7, |2 3. |2 V + 15(? 



f2> 



dt V 2 ' 2 
The RGE for the scalar-potential couplings are 

2 



1/4, 



2/5- 



(25) 
(26) 
(27) 



16n 



lQn 2 



lQn 2 



lQn 2 



lQn 2 



16tt 2 



dt 



dA2 
dt 



dAa 
dt 



dA 4 
~dt 



dAs 
dt 



dAe 
dt 



2A\( + + (A 4 + MY + K + (A 5 + A 8 ) z + 4 



+ (4|y 3 | 2 -C) Ai 



9<? 4 , 3g 2 g' 2 , 3 5 



,/4 



+ 



-2 



24A 2 + A 2 + (A 4 + A 7 ) 2 + A" + (A 6 + A 9 ) 2 + 4 

/ , ,2 ^\ , 9# 4 3#y 2 3#' 4 A 
+ (8 |y 4 | 2 - C) A 2 + -|- + ^f- + -f- - 4 

24A 2 + A 2 + (A 5 + A 8 ) 2 + A 2 + (A 6 + A 9 ) 2 + 4 



+ (8|ys| 2 -C) A 3 + ^- + 



9<? 4 , 3g 2 g> 2 , 3<?' 



/4 

+ — 4 



Aio| 2 + 4|A n | 2 



2/s I , 

Aio| 2 + 4|A 12 | 2 
2/4| 4 , 

An| 2 + 4|A 12 | 2 
2/5 1 4 , 



8 ' 4 

(Ai + A 2 ) (12A 4 + 4A 7 ) + 4A 2 + 2A 2 + 4A 5 A 6 + 2 (A 5 A 9 + A 6 A 8 ) - 8 : \ 

f 



(28) 
(29) 
(30) 

|2 
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+ ( 2 N 2 + 4M 2 -c)a 4 + ^-^ + 3 ^ 4 



(31) 



2 ' 4 ' 

(Ai + A 3 ) (12A 5 + 4A 8 ) + 4A 2 + 2A 2 + 4A 4 A 6 + 2 (A 4 A 9 + A 6 A 7 ) + 8 |A n | 2 



+ (2 



| 2 + 4|, 5 | 2 -C)A 5 + ^-^ + ^, 



(32) 



(A 2 + A 3 ) (12A 6 + 4A 9 ) + 4A 2 + 2A 2 + 4A 4 A 5 + 2 (A 4 A 8 + A 5 A 7 ) + 8 |A 12 | 2 
+ (4 M 2 + 4 |y 5 | 2 - C) A 6 + % - ^f- + - 8 |^ 5 | 2 , (33) 



16tt 



2 dA? 
dt 



(4Ai + 4A 2 + 8A 4 + 4A 7 + 2 |y 3 | 2 + 4 \y 4 \ 2 - C) A 7 

7 



16tt- 



16tT 



16tt' 



16tH 



> dAs 
dt 

j dAg 

dt 

dAio 
dt 

dAn 
dt 



+2A 8 A 9 + 32|A 10 | 2 + 3( ? V 2 , 

(4Ai + 4A 3 + 8A 5 + 4A 8 + 2 \y 3 \ 2 + 4 |y 5 | 2 - C) A 8 



+2A 7 A 9 + 32|A 11 | 2 + 3( ? V 2 , 



16tt 



2 dAi2 
dt 



where 



(4A 2 + 4A 3 + 8A 6 + 4A 9 + 4 \y A \ 2 + 4 \y 5 \ 2 - c) A 9 

+2A 7 A 8 + 32 |A 12 | 2 + 3g 2 g' 2 - 8 \y 4 y 5 \ 2 , 

(4Ai + 4A 2 + 8A 4 + 12A 7 + 2 |y 3 | 2 + 4 |y 4 | 2 - C) A 10 

+4AhA; 2 , 

(4Ai + 4A 3 + 8A 5 + 12A 8 + 2 \y 3 \ 2 + 4 \y b \ 2 - C) A n 
+4A10A12, 

(4A 2 + 4A 3 + 8A 6 + 12A 9 + 4 \y A \ 2 + 4 |y 5 | 2 - c) A 12 
+4At An - 4yf y 5 2 , 

_/2 



(34) 

(35) 

(36) 
(37) 
(38) 

(39) 
(40) 



C:=9g 2 + 3g' 

The reason why no fourth-order terms in the Yukawa couplings and in the gauge couplings 
appear in the RGE for A10 and An is that the condition A10 = An = may be enforced 
through an additional U(l) symmetry: <f>i — > e l0t <f>i, — > e~ ta e^, where a G R. 

We next write down the RGE for the coupling matrices of the effective neutrino mass 
operators. They are 



16^ 



16tt' 



16tt 2 



16tT 



dt 

d^2) 
dt 

d^3) 
dt 

dt 

^(32) 
dt 



(-3g 2 + 4A X + 2 |y 3 | 2 ) «<"> + 4A* ^^ 22 ) + 4AJ 1 k< 33 > 
+ {/<( 11 ),P-2P 1 }, 

(-3/ + 4A 2 + 4 |y 4 | 2 ) ^ + 4A 10 k( 11 ) + 4A* 12 k< 33 > 
+ {k^ 22 \ P - 2P 2 } - 2 (k^P 23 + iW 32 *) , 

(-3s 2 + 4A 3 + 4 |y 5 | 2 ) K ( 33 ) + 4A n « (11) + 4A 12 *; (22) 
+ { K ™, P - 2P 3 } - 2 (^ 32 )p 32 + P 32 ^ (23) ) , 

(-3/ + 2A 6 + 2 M 2 + 2 \y,\ 2 ) ^ + 2A 9 k< 32 > + {k™, p} 



(41) 



(42) 



(43) 



-4k (22) P 32 + 2P 32 k {22) - 4P 23 /t {33) + 2k (33) P : 



23 



+ 2 



^ 23 ),P 2 -P 3 ] -2(^ 32 )P 3 + P 2 ^ 32 )), 

(-3^ 2 + 2A 6 + 2 |y 4 | 2 + 2 |y 5 | 2 ) K ( 32 ) + 2A 9 k< 23 > + {«<»>, P} 

-Ak^P 23 + 2P 23 /t (33) - 4P 32 «( 22 > + 2k^P 32 
+2 [k™, P 3 - P 2 ] - 2 ( K ( 23 )P 2 + P 3 k^) , 



(44) 



(45) 



where {R, S} and [R, S] denote the anticommutator and the commutator, respectively, of 
the matrices R and S. Moreover, we have defined 



Pi 


:= diag 


(M 2 , 0, 0), 


(46) 


Pi 


:= diag 


(0, IvaI 2 , M 2 ) , 


(47) 


Ps 


:= diag 


(0, \y 5 \ 2 , N 2 ) , 


(48) 


P23 


:= diag 




(49) 


P32 


:= diag 


(0, 2/42/5, -vIvb)- 


(50) 



Notice that the matrix P in @ is equal to (P 1 + P 2 + P 3 ) /2. 



4 Predictions of the Z2 and D 4 models 

The Lagrangian of neutrino Majorana masses is 

1 



a,f3=e,fi,T 



•^Majorana 

where M. v = M.^. Taking b = d = 2 in (j^J), it is clear that 



(51) 



(52) 



Therefore, if we denote the VEV of 0° by Vi, then the neutrino Majorana mass matrix 
M. v is given by 



1 



^3 



i=l 



+ K< 32 > 



(53) 



since = = = = in the Z 2 and D 4 models. This is valid at all scales 
t. 

In general one may write [2U| 

/ X A(l + e) A(l-e)\ 
M u = A(l + e) B(l + e') C . (54) 
U(l-e) C B(l-e')J 

We already know that the form of the flavour coupling matrices is described by 
while /«( 23 ) = z^ 32 ) is decribed by ([24)1 . Therefore, 



2 ,,(«•) 



while 



i=l 



:= K ( 23 ) + «< 32 > 



.4 








C 


1 


C 









Cl 


-Cl 


C\ 


c 2 





-Cl 





-c 2 



(55) 



(56) 



with 



eA = 2v 2 v 3 ci, 
e'B = 2v 2 v 3 c 2 . 



(57) 
(58) 



Once again, all this is valid at any scale t. 

In both the 7j 2 and D4 models, the symmetry Z 2 aux ^ inverts the signs of the right- 
handed-neutrino fields which are present above the seesaw scale. Hence, only the doublet 
<pi has Yukawa couplings to those fields, above the high scale. This implies that M. v (to), 
where to := lnm#, originates solely from the VEV of cf>\ |1H I12|. Therefore, at the seesaw 
scale [Hj 

K^(t ) = M u (to)/(2vf), 
«(«) (t ) = for all other (ij). 



(59) 



We conclude that M. v (t 

\T ■ 



has the same form as i.e. M. v (to) is of the form (|23p. 

Clearly then, (0, 1, — 1) J is an eigenvector of M. v (t ) and therefore, at the seesaw scale, 
the predictions (fT3)l hold. 

At any other scale, though, the matrix in (156)) is not zero. Thus, for any t < t , 
M. v (t) is not fi <-» r-symmetric. This fact renders the predictions (J 14)) inexact for any 
scale other than the seesaw scale. In j2D] it has been shown that, if one assumes the 
parameters e and e' in ()54|) to be small, then, to first order in those parameters, one has, 
instead of (JHJ), 



V, 



c3 



~2 I es l2 m 2 + e s l2 m 3 — e m 2 — e m 3 
rn 3 — rn 2 



+ 



S12C12 /_ 2 



2 ( ec 12 m 1 + e c 12 m 3 + em 1 + e m 3 



m 3 — mf v 



(60) 



cos 20. 



23 



Re 



2c 2 



12 



m, — nv 



( es i2 - e ') (rh 2 + m 3 )* 



2s 2 



( CC 12 + e ') ("H + m 3) 



m 3 — mf v 



where 



[m 1 -m 2 )e, 
rhis\ 2 + rh 2 c\ 2 + m 3 , 
2 6 ' 



and 



mi 
m 2 



m ie - 2 ^, 



-2lCT 



(61) 

(62) 
(63) 



(64) 
(65) 



Here, m 1; m 2 , and m 3 are the (real, non- negative) neutrino masses, while p and a are 
Majorana phases. In (|6l} )l -fH3 ^ . S\ 2 and Ci 2 are the sine and cosine, respectively, of the 
solar mixing angle. We are using the standard parametrization EI! for the lepton 
mixing matrix U; that parametrization fixes the significance of the phase of (|60j) . 
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We see in (|6U|) and (|61j) that the deviation from the predictions (J 14)) is numerically of 
the same order as the (small) parameters e and e', except in the case of quasi-degenerate 
neutrinos with common mass mo; in that case, an enhancement factor mg/Am^ appears, 
where Am 2 tm = |mj — m 2 | « |m| — m 2 | is the atmospheric mass-squared difference. 



5 An approximate solution of the RGE 

In this section we assume the Yukawa couplings y 3 , y 4 , and y 5 to be very small. This 
assumption is reasonable when one considers the Z 2 or D4 models with their minimal 
content of three Higgs doublets, since in that case J2i=i \ v i\ 2 must be equal to (174 GeV) 2 , 
and therefore all the \vi\ are in principle much larger than the charged-lepton masses. We 
estimate, to first order in the Yukawa couplings squared, the deviation of M. v (ti), where 
t\ := In ruz, from the \l <-> r-symmetric form. We thus find out the likely magnitudes of 
U e s and cos 2# 23 at the electroweak scale. 

The starting point is (|4l" |) -(|4l) j) . We drop all the Yukawa couplings from (JUjHfjEI) ; 
in f)44j) and (|45)1 . on the other hand, we keep the terms with P 23 and P 32 , since it is those 
terms which induce corrections to the fi <-> r-symmetric form of M. v (to)- We thus obtain 



16^ 



16vr" 



16vr 



dt 
d^ 22 > 

dt 
d^ 33 > 

dt 



and 



16vr 



dn^ 



dt 



3g 2 + 4Ai) + 4A* «: (22) + 4A* 1 k:( 33 >, 
3g 2 + 4A 2 ) + 4A 10 k (11) + 4A^ 33 \ 
(-3/ + 4A 3 ) k (33) + 4A llK ( n ) + 4A 12 k( 22 ), 

V + 2A 6 + 2A 9 ) ^ - 2 {P 32 , «: (22) } - 2 {P 23 , ^ (33) } , 



where 



K 



(23) 



+ k^ 32 ) as in (JSEJ). Using (gSJ) and (jSUJ), we obtain from fU) that 



167T' 



dci 

~d7 

dc 2 



(~3g 2 + 2A 6 + 2A 9 ) ci - 2 [^i/s/q? + 2/42/^ 



(33) 
12 



16vr^ ^ = (-3# 2 + 2A 6 + 2A 9 J c 2 - 4 [y&sKja + ViVl^n 



-(33) 



(66) 
(67) 
(68) 

(69) 

(70) 
(71) 



Formally we can write the solution of the coupled differential equations (|66j) - (j68|) as 



«W (t)=£T(M')^ fe) (f). 
3=1 



(72) 



All the matrix elements of evolve according to the same operator T(t,t'). Since at 
the scale t only is non-vanishing, 

^(t)=T(t,t ) il K^(t ). (73) 

Formally, the solutions to (J7(Jj) and (J71j) are easily written down. Defining 



S c (t) : = exp 



16tt 2 



^dt'(-3(7 2 + 2A 6 + 2A 9 ) (0 



(74) 
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and taking into account that c\ (t ) = c 2 (to) — 0, one has 

ci(t) = ~S e (t) fdt'S^^lyty^+y^K^it 1 ) 

07T Z Jt 1 J 

K (11) (tn) ^ 

' J S c (t)/ dt'S; l (t>)[v(t>)T(t>,t ) 21 + v*(t')T(t',t ) 31 ], (75) 



8tt 2 

K 22 ^ (to) 



to 



'to 

— . ~t i / i a i u . 'i 

4tt 2 



S c (t) / dt' S" 1 (0 [i/ (0 T (f, t ) 21 + (f ) T (f , t ) 31 ] , (76) 



to 



where we have used (J7HJ) and defined v := y\y^. 
We now make use of (J55J1 and (JSTj) to write 



V2V3C1 (ti 



V2V3C1 (ti) 1 

Eti^(*i,*o)u 4^ (t )' 



(77) 



Similarly, 



, , s_ V2V3C2 (h) I , 

[1) ZUvfT^Ma^itoY [ 



Putting (|75|) - (|75|) together, we conclude that 



V2V3 



il 



e (*0 = - c 9V ^3 ^ ~ ^ (*0 / " d*' S; 1 (f) [1/ (f ) T (f, t ) 21 + u* (f) T (t\ t ) 31 ] 

*0 



(79) 

while e ; (ti) = 2e (ti). From (JoTlj) and ()61|) one can derive that 



2e 



/mj + m 3 m*+m 3 \ 

U e3 = 2m 3 ci 2 si2 —5 n H o 2 Ree > 

\ mo — mi mi — mi I 

2 2x (8°) 

cos 2^3 = 2 sl J A \ +m i +cl J A \ +m i M 
V mf — m 3 TO2 — m 3 / 



Now we want to estimate the maximum possible order of magnitude of e (ti) by us- 
ing (|79|). The length of the integration interval of t' is t — ti = In {m^jmz) ~ 10. The 
functions S c (t) and S 1 ^ 1 (t') are of order 1 since, in (|74p. A 6 / (167T 2 ) and A 9 / (167T 2 ) are 
necessarily small. The functions T(t',t ) 21 and T(t',t ) 31 are governed by A10 and An; 
in any case, T (t', to) 21 /T (ti, to) a and T (t', t ) 31 /T (ti, to)ji should be ^ 1. Similarly, 
V2V3/V? should not be larger than 1. We conclude that 



| e( tO|~^l~l^. (si) 
Svr 2 10 



Even if we allow for rather small VEVs, 1 2/4,5 1 cannot be larger than 0.1. We thus have 
the generous upper bound |e (ti)| < 10 -3 . 
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Equation (|8U|) tells us that the only chance to have a non-negligible is in the case 
of a degenerate neutrino spectrum. Let us consider the extreme case of a common mass 
m = 0.3 eV |2lj. Since Am^ tm ~ 2xlO~ 3 eV 2 , in that case we have |[/ e3 | ~ 100 |Ree| <0.1. 
Here we have used Majorana phases p = o = for simplicity. That choice of mo is 
indeed extreme; if take m = 0.1 eV instead, then the upper bound becomes one order 
of magnitude smaller, due to the quadratic dependence of \U e s\ on m in the case of a 
degenerate neutrino spectrum. In any case, we expect \U e3 \ and | cos 2^1 to be no larger 
than 0.1 in our model, but most likely they are two or more orders of magnitude smaller. 

6 Summary 

In this paper we have computed the RGE for the dimension-five neutrino mass operators 
in the multi-Higgs-doublet SM. Thus, the main result of this paper is (JJJ), which describes 
the evolution of the coupling matrices of the mass operators in the SM with an arbitrary 
number of Higgs doublets. We have argued in favour of the usefulness of (j2J) by citing 
models for lepton mixing which have been constructed in the framework of the multi- 
Higgs-doublet SM. 

As an application of our RGE we have considered the Z 2 model of [ITj and the -D4 
model of which — from the field-theoretical point of view — are identical below the 
seesaw scale. The predictions (fT3j) of those models hold at the seesaw scale and we have 
used the RGE to estimate the corrections to those predictions which appear due to the 
evolution of the coupling matrices of the dimension-five neutrino mass operators down to 
the electroweak scale. We have found that those corrections are in general negligible, with 
the possible exception of a degenerate neutrino mass spectrum with a rather large common 
mass m > 0.2 eV. In that case, s 2 3 could be as large as 0.01 and be within the sensitivity 
of the planned long-baseline neutrino experiments jU]. On the other hand, even in the 
degenerate neutrino case the deviation of #23 from 7r/4 will be hard to uncover in those 
experiments [22], since they will be sensitive to the parameter sin 2 2# 23 = 1 — cos 2 26* 23 and 
we have estimated cos 2 2# 23 ~ 0.01 in our models. Thus, with respect to the experiments 
presently envisaged the models discussed here have the following properties: should a 
non-zero s 2 3 be discovered, then the neutrino mass spectrum must be degenerate; while 
deviations from sin 2 2# 23 = 1 should be invisible. 

A Vertex corrections to the neutrino mass operator 

The last two lines of (J7|) originate in vertex corrections of the type displayed in figure 
In this appendix we show how we have arrived at those two lines. Perturbation theory 
yields the expression 




x 




k=l e=l 



1=1 f=l 



(Al) 
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where d is the dimension of space-time and x, xi, are space-time points. We have left 
out the flavour indices in (jAlj) . The gauge-S77(2) indices are a,b, . . . , f. The symbol T 
denotes time ordering. When computing (jAlj) . the field £r must be contracted with 
As for D-Le-, it may be contracted either with D^ a or with D^ c ; it is easy to see that both 
possibilities yield the same contribution to the RGE of k^' — this fact explains the factor 
2 in the second line of (j7j) . In the following we compute explicitly the case where Dl s is 
contracted with D Lc . 

For the contraction of <p[f there are also two possibilities: one may contract it either 
with (fimb or with 4>nd- We consider the second possibility first. We use dimensional 
regularization with minimal subtraction. In the evaluation of (jAlj) we only need the pole 
terms in e = 4 — d. The computation is straightforward and we arrive at 

i n H 2 

2~ E E Dl a C- 1 [K^Y^Y n }D Ld e ab ^ mb e d ^ kc , (A2) 

k,m,n=l a,b,c,d=l 

where all the fields are now meant to be at the same space-time point. The minus sign 
from the factor (— i) 2 in (jAlj) has been removed through the interchange of the indices c 
and d in e cd . If, instead, we contract <f>h with <p mb , then we get 

i n H 2 

E E D^C-^K^YiYJl^^M^. (A3) 
lo7r z e , , , , 1 L J 

k,m,n=l a,b,c,a=l 

The SU(2) structure of (jA3j) is different from the one in (j4*jl . Therefore, we need to apply 
the identity 

£ ab £ cd + £ ac £ db + £ ad £ bc = g ( A4 ) 



to (jA3j) . obtaining 

1 n H 2 

E E i^" 1 [« (mn) n f ^] (A5) 

^* ^ fc,m,n=l a,6,c,ci=l 

The terms in the last two lines of ((7j) are obtained from (jA2j) and (jA5j) as follows. Firstly, 
one substitutes the factor 1/e by — 1 [7j. Secondly, the indices k and m — in (jA2j) — or k 
and n — in (jA5j) — of the scalar doublets must be replaced by % and j; the contribution 
to the beta function of is then given by the flavour matrix in between the lepton 
doublets. If, in each expression, the first Higgs doublet is labeled i and the second one is 
labeled j, then one obtains the terms in the second line of (J7J): — Yin ^ m ^YjY n from (jA2j) . 
Ym ^ m ^Y^Y m — Ym K^YjYmj from (jA5j) . Reversing the role of i and j leads to the 
terms in the last line of (J7|). 
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